The besj1() function returns the j1st Bessel function of its argument, i.e it returns the first order Bessel function of the first kind. besj1() expects its argument to be in radians.